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I review some important aspects of the structural and statistical properties of the nearby 
£C) . X-ray galaxy cluster population, discussing the new constraints on mass profiles, the mass- 

if temperature relation, and the entropy of the intracluster medium which have become available 

ON from recent X-ray observations. 

On 

o , 

O 1 Introduction 

o 

Massive systems of galaxies, gas and dark matter, located at the intersection of tenuous, large- 
scale filaments, galaxy clusters are the nodes of the highly structured Universe we see today. 
Formed from the collapse of initial density fluctuations, clusters are built hierarchically and 
constitute an evolving population. The evolution of structure in the Universe - and thus that 
of clusters - is strongly dependent on the cosmology, so that the statistical properties of the 
cluster population are a test both of cosmology and of structure formation theory itself. 

While clusters were initially discovered as over densities of galaxies on optical plates, our 
present understanding of their composition suggests that the galaxies themselves constitute 
only ~ 5 per cent of the total mass, with most of the mass (~ 80 per cent) in the dark matter 
component. The dominant baryonic component is the hot, X-ray emitting intracluster medium 
(ICM), a rarefied plasma at several million K, reflecting the great depth of the potential well 
in which it lies. We are currently in the privileged position of having several world-class X- 
ray observatories with which to examine the structural and scaling properties of the ICM in 
galaxy clusters. The arcsecond spatial resolution of Chandra is complemented by the grasp 
of XMM-Newton, allowing spatially-resolved observations to be performed with unprecedented 
precision 3 . 

In this paper I begin with a brief review of the expectations for structure and scaling from 
the starting point of simple gravitational structure formation. I then compare these simple 
expectations with recent observations of nearby clusters from XMM-Newton and Chandra. I 
first show how X-ray observations can be used to place strong constraints on the dark matter 
component, and then, passing via the M-T relation, I take a look at how measurements of the 
entropy have cast light on the effect of non-gravitational processes, such as galaxy formation 
and cooling, on the X-ray properties of the cluster population. 



a In addition, the recently- launched Suzaku X-ray satellite has a much lower instrumental background, holding 
the promise of spatially resolved observations to larger cluster-centric distances. 



2 Expectations for structure and scaling 



It is useful to have a baseline prediction for the X-ray properties of the cluster population. Cluster 
formation is driven by the gravitational collapse of the dominant dark matter component. The 
ICM collapses with the dark matter, so that, to a first approximation, the gas properties should 
be determined entirely by the processes involved in gravitational collapse, i.e ., shock heating and 



compression. The baseline model rests on several simple assumptions (e.g.. lEke et al 



The internal structures of clusters of different mass are similar. 

All clusters identified at a given redshift correspond to a given characteristic density, and 
the characteristic density scales with the critical density of the Universe, i.e., 

^ = <p> = ^W*);* -zoo. (i) 

^200 6 

The ICM evolves in the gravitational potential of the dark matter; the gas mass fraction 
is thus constant, 



f - M e as (2) 

Time-averaged, the ICM is in approximate hydrostatic equilibrium, allowing application 
of the Virial Theorem, giving: 
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where T is the mean cluster X-ray temperature and fir is a normalisation constant depen- 
dent on internal structure. 

The above equations imply that power-law relations exist between various X-ray properties, Q, 
and the mass M (or a mass proxy) and the redshift, z, such that Q cx A(z)M a . The evolution 
factor, A(z), is due to the evolution of the mean density, which varies with the critical density, 
p c (z) cx h 2 (z). For example, the gas and total mass scale as M gas cx M 200 oc br x {z)i*l 2 , 
and, assuming bremsstrahlung emission, the X-ray luminosity scales as Lx oc h(z)T 2 . The 
assumption of structural similarity further implies that, once scaled appropriately, the radial 
profiles of various cluster quantities (e.g., gas density p, entropy S, temperature T, total mass 
M) are similar. 

Clusters created in hydrodynamical simulations of purely gravitational structure formation 
are structurally similar and follow these scaling relations (e.g.. B rvan fc Norm ail Il99?i ). showing 



that the assumptions underlying this simple spherical baseline model are applicable in a hierar- 
chical cosmological context?. In comparison with observations, the study of both the correlations 
between global properties (the scaling relations) and the structural properties (radially-averaged 
profiles of various quantities) is essential for the best understanding of cluster formation and 
evolution. 



''Note that there is a certain amount of intrinsic scatter in the relations, which is related to cluster dynamics. 
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Figure 1: Left: In tegrated mass profiles of 1 clusters in the temperature range 2-9 keV, observed with 
XMM-Newton (from Ipointecouteau et aHl2005l) . For each cluster, the mass is scaled to M200 and the radius to 
i?200- The black line corresponds to the mean scaled best-fitting NFW model;, the dashed lines the associated 
standard deviation. Right: Concentration parameter vs mass relation for NFW model fits to the total density 
profiles of 12 clusters in the temperature range 0.7 - 9 keV, observed with Chandra ffrom fvikhlinin et all 2006). 
The solid line shows the average concentration of CDM haloes from simulations; the dotted lines show the 

associated 2a scatter. 



3 Dark matter constraints 



3. 1 Mass profiles 



An inescapable prediction from more than a decade's worth of numerical simulations is the 
existence in haloes of a semi-universal cusp ed dark matter dens ity profile. The most well- 
known variant is the Navarro- Prenk- White ( Navarro et al. . Il997l . NFW) profile, Pdm{ t ) oc 
[(r/r s )(l + r/r s ) 2 ] -1 , where the concentration parameter, eg = Rs/r s , is a measure of the 
halo concentration. The profile is not strictly universal since there is a weak dependence of C5 
on mass, Ma , due to the fa c t that smaller mas s haloe s generally form earlier, when the Universe 
was denser (Bul lock et all . 12001 : iDolag et all l2004h . Although the total mass density profile 
can be measured with both X-ray and optical observations, in the latter case, stacking is often 
required to overcome the paucity of data on individual objects if qu antitative measurements of 
the mass profile are to be obtained (e.g., Biviano Girardl 2003). in contrast, while X-ray 
observations allow measurements of individual mass density profiles, the method relies on the 
assumption that the gas is in hydrostatic equilibrium. 

X-ray observational confirmation of the exist ence of a cusped p r ofile o n cluster scales started 



to become available with ASCA observations (Ma rkevitch et al. although the spatial 



resolution of these observations, in particular of the temperature profiles, was not optimal. 
More recent observations with XMM-Newton and Chandra have allowed t he mass profiles o f 



relaxed clusters to be probed from the ver y core regions (~ O.OlijWi. e.g.. iLewis et al 



right out to 0.7 - 0.8i??nn (i.e., ~ Rwn; IPratt fc Arnaudl 120021: Ipointecouteau et al 



2003) 



2004; 



Vikhlinin et~all l200fil : IVoigt fc Fabianl 200g|) . All of these observations point to the existence of 
a quasi- universal cusped profile as predicted from CDM simulations . There is no evidence for 
a core in even the highest-resolution observations ((Lewis et al ., 2003), effectively ruling out the 
possibility of self-interacting dark matter on cluster scales. 



3.2 Concentration parameters 



Observations of small samples of clusters and gal axy groups have also al l owed the universality 
of mass and total density profiles to be prob ed jPratt fc Arnaudl 1200.4 iPointecouteau et al 



Vikhlinin et all 2006; Humphrey et all l200fih . The left-hand panel of Fig. ^ shows the 
integrated mass profiles of 10 clusters in the temperature range 2-9 keV, scaled by the virial 
mass M200 and radius i?200- The mass profiles can all be adequately fitted with NFW models 
and are clearly similar once scaled in this manner. The right-hand panel of Fig. ^ shows the 
measured concentration para meters from fits t o the total density profile of 12 clusters in the 



of Dolag et al 



te mperature range 0. 7-9 keV ((Vikhlinin et all |2Q06), compared to the theoretical predictions 



(2004). The observed concentration parameters are in good agreement, both 
in absolu te value and dispersion, wi th the theoretical predictions. Similar results have been 



found bv IPointecouteau et al.l (J2005) for another sample of clusters, and bv iHumphrev et al 
(|200fih for a sample of early type galaxies and galaxy groups. In the latter case, it was found 
that inclusion of a stellar mass component was necessary to give acceptable fits to the data, 
since this component starts to dominate the mass budget in the very inner regions (see also 

Such results are encouraging and have allowed constraints to be put on 



the nature of the dark matter itself on cluster scales, suggesting that our understanding of the 
dark matter collapse is correct down to the scale of individual giant galaxies. 

Finally, it should be pointed out that all of the above results were obtained only for mor- 
phologically relaxed systems. Attempts to model the to tal mass dens i ty profiles of less relaxe d 
systems with the NFW profile have resulted in poor fits feelsole et all 1200.4 IPratt et all 1200.4 . 
This suggests that in these cases either (i) the gas in not in hydrostatic equilibrium, or (ii) that 
the underlying DM density profile is not well described by an NFW model (the original NFW 
model was proposed specifically for systems with a relatively high degree of relaxation). 



4 The M—T relation 



Models of structure formation predict the space density, distribution and physical properties of 
the cluster population as a function of mass and redshift, so that, from the theoretical point of 
view, the mass of a cluster is its most fundamental property. However, X-ray cluster surveys 
yield the cluster space density and distribution only in terms of X-ray observables such as the 
luminosity, Lx, or the temperature T. Scaling relations linking the mass to the X-ray observables 
allow the full use of X-ray cluster survey information in the statistical sense. Knowledge of mass 
observable relations and the associated scatter about these relations is thus crucial in the use of 
clusters to probe cosmology. 

Assuming hydrostatic equilibrium and spherical symmetry, the mass of a cluster can be 
recovered from spatially-resolved X-ray observations of the density and temperature profiles. 
The M—T relation is one of the fundamental relations for linking the observed gas properties 
with theory, since the temperature is expected to be closely related to the mass via the Virial 
theorem. It is also of topical interest, since the normalisation of this relation has a direct bear- 
ing on the measurement of the present-day value of ag using cluster abundance measurements, 
and, until the WMAP 3-year data, cluster measurements of erg were thought to be unusually 
low. Samples from ROSAT and ASCA, though large, either assumed that clu sters were isother- 
mal (|Horner et all Il999l : IXu et all l200ll : ICastillo-Morales & SchindlerL l2003h or had relatively 
poor temperature profile reso l ution and /or relied strongly on extrapolation t o derive the virial 
mass ( Nevalainen et al. . 2000; Fi noguenov et al. . 200ll : Sanderson et al. . 2003). As a result, his- 
torically, there has been little consensus on either the slope or the normalisation of the M—T 
relation. 

The situation is much improved with recent results from XMM-Newton and Chandra. The 
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Figure 2: Left: The X- ray mass-temperature relation from obse rvations obtained with 

XMM-Newton iArnaud et alJ . 120051 circles) and Chandra iVikhlinin et all I2006L stars). Some clusters 
were observed with both satellites, and are plotted twice. The long-dashed line is the best-fitting relation to 
the XMM-Newton data points, the dashed line is the best fitting relation to the Chandra points, and the solid 
lines is the best fitting relation to all data points. See text for details. Right: Comparison with results from 
several numerical simulations. The inclusion of more realistic physical processes in more recent simulations has 

improved the agreement between observation and theory. 



left-hand panel of Figure [21 shows t he mass-temperature data from two cluster samples observe d 
independently with XMM-Newton ( Arnaud et all l2005h and Chandra (jVikhlinin et al.l .Eo06y. 
The various lines show the results of fitting a power-law function of the form h(z)M = Ax 
(kT/5 keV) a to different data sets. The fits to the individual XMM-Newton and Chandra data 
sets give results which are entirely consistent within the la errors for both slope and normali- 
sation. A BCES fit to the total data set yields the relation 



log[/i(z)M 5 oo/lO 14 M ] = (0.57 ±0.02) + (1.59 ± 0.08) log [kT/5 keV]. 



(4) 



The slope is thus very close to the expectations from the simple scaling analysis in Sect.[2]above. 

A second historical point of contention is the disagreement between the observed M-T nor- 
malisation and that derived from numerical simulations. Now that there is good agreement 
between measures of the observed normalisation, the time is ripe for a new comparison with 
simulations. The right-hand panel of Figure 121 shows the data an d best fitting power law de- 
scribed above, compared to the relations found in the simulations of Evrard et al. dl996h . which 



onl y include gravita tional processes, and the more recent simulations of hwar,^ B 



a Aravtsovetakl »), which include additional physical processes such as radiative cc^iuTg 
and feedback from supernovae. The disagreement with the gravitation-only simulations is well 
documented and for the present data set amounts to a difference in normalisation of 40 per cent 
at 5 keV. The agreement has clearly improved with the introduction of more realistic physical 
processes into the simulations. 

The remaining discrepancy is likely due to a combination of two effects. Firstly, it is unclear 
how simplifying assumptions of the X-ray analysis have an impact on the final mass determi- 



"Similar results were found at -R2500 for a smaller sample of 6 systems bv lAllen et alJ i200 



nation. In particular, the validity of the assumption of hydrostatic equilibrium has often been 
called into question. Observation s of Coma suggest a lo wer limit of ~ 10 per cent of the total 



ICM pressure in turbulent form ( Schuecker et al. . 20041 ). However, Coma is a merging cluster 



and at present no observational constraints exist for the amount of non-thermal pressure support 
in clusters which observers would call 'relaxed'. Secondly, t he value of the M—T norma l isatio n 
can vary by up to 50 per cent depending on the simulation ( Arnaud et al. , 200,4 Henrv , 2004h , 
and the normalisation has been creeping downwards over time as more complex physics has 
been incorporated into the simulations. Current simulations are still unable to reproduce not 
only the observed scaling relations (e.g., the Lx~T relation), but also the structural properties 
of clusters (e.g., temperature profiles in the central regions), clearly pointing to the possibility 
that certain physical processes may be missing. Below I discuss how X-ray observations can give 
an idea of what other physical processes may be in play. 

Finally, it should be noted that the above results were obtained for morphologically relaxed 
cluster samples only. Since the precision on the estimation of cosmologic al parameter estimation 
depends on the precision to which the scaling relations are known (e.g.. iHenrv . l2004h . accurate 
estimates of the dispersion around the best fitting relations are needed for the whole cluster 
population. This is true even for the case of self-calibration of extremely large samples. The 
challenge now is to derive accurate mass-observable scaling relations for representative cluster 
samples. Since the X-ray method is clearly invalid for clusters which are not in hydrostatic 
equilibrium, this will require inter-calibration of the available mass estimators (X-ray, optical, 
lensing), preferably using both observed and simulated data sets. 



5 Gas physics 



X-ray data have been telling us for nearly two decades that the simple scaling relations outlined 
in Sect.|2]do not describe the observed cluster population. The most well-studie d example is that 
of the Lx-T relation, which has been found by many studies to be Ly oc T 3 ( Edge fc Stewartl . 
19911 : lArnaud fc EvrardL Il999l : iMarkevitcrj . ll99Sl : I Osmond Sz PonmanL l2004h , corresponding to 
a suppression of X-ray luminosity in poor systems relative to simple expectations. This suggests 
that physical processes other than gravity alone are affecting the properties of the ICM, and 
that poorer (i.e., lower temperature) systems are proportionally more affected. These non- 
gravitational processes are likely linked to radiative cooling and to heating processes associated 
with galaxy formation. 

The entropy affords the most direct way of investigating the properties of the ICM gas. The 
entropj^ is generated in shocks as the gas falls into the deep potential of the cluster halo. Entropy 
is conserved in any adiabatic rearrangement of the gas, and the gas will always rearrange itself 
so that the entropy increases outwards. This property is illustrated in the left-hand panel of 
Fig. 13 which shows the entropy profiles of 40 clusters spa nning a mass ran ge of a factor of ten, 



from simulations including only gravitational processes (|Voit et al 
been scaled by the characteristic entropy of the halo 



2005). The profiles have 
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(5) 



where f b is the baryon fraction, and the radii have been scaled to the virial radius i?200- It can 
be seen that the profiles all coincide, and that, outside the central regions (where gravitational 
softening becomes important in the simulations), a power-law model with K oc R is an 
adequate description of the trend with radius. In the right-hand panel of Fig. |21 I show the 



d Here I follow the established convention and refer to the quantity S = K — kTn e 2 ^ 3 as the entropy. This is 
related to the true thermodynamic entropy via s = lnK 3 ^ 2 + constant. 
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Figure 3: Left: Entropy profiles of 40 clu sters spann i ng th e mass range 0.2 — 8 x 10 h Mq from simulations in- 
cluding only gravitational processes (from IVoit et all l2005h . The profiles have been scaled using the characteristic 
entropy of the halo K200 and the virial radius -R200 (see text for details). Right: Entropy profiles of 10 clusters 
spanning th e mass range 1 - 12 x 10 14 /ifo M @ (~ 2 [blue] - 8keV [red]), observed with XMM-Newton (from 
iPratt et , all I2OO6I . The profiles have been scaled in exactly the same manner. The solid line is the best-fitting 
power law description of the radial dependence of the entropy for the simulated clusters. 



entropy profiles of ten morphologically relaxed clusters spanning th e temp erature range 2-10 
keV (~ 10 14 - 1O 15 M ), observed with XMM-Newton (|Pratt et all l200fl L The clusters are 



colour-coded so that poor (cool, low mass) clusters are blue and rich (hot, high mass) clusters 
are red. They have been scaled in exactly the same manner as the simulated profiles, assuming 
fb = 0.14 (Ofe/1 2 = 0.02 and Q m = 0.3), and the solid line denotes the best-fitting power-law 
entropy-radius relation to the scaled simulated clusters. It can be seen that, while the profiles of 
rich clusters are in good agreement (both in slope and in normalisation) with the prediction from 
gravitational entropy generation, poor clusters have a systematically higher entropy throughout 
the ICM. 

We can get an idea of the dependence of the offset with temperature (or mass) by looking 
at the entropy measured at a certain fraction of the virial radius. The entropy-temperature 
S-T (or K-T) relation was extensively investigated using data from the previous generation 
of X-ray satellites. In the simple baseline model outlined in Sect. [21 the gas mass fraction 
is constant, implying a constant gas density and thus a simple linear scaling of the entropy 
with temperature S oc h(z)~ 4: / 3 T ' , or m ass S oc /i(z) _2 / 3 M 2 / 3 . Initial measurements of the en- 
tropy at 0.1 i?200 ( Ponman et allll999l ) showed a levelling-off of the S-T relation towards lower 



temperatures, which was subsequently interpreted as a limiting value to the central entropy 
( Llovd-Davies et all EoOoT). However, dee per investigation with larger samples observed with 



ROSAT and ASCA (Ponma n et all I200.T) . and smaller samples observed with the newest gen- 
eration of satellites ((Pratt et all 120061 ) have shown that there is in fact a continuous power-law 
relationship between the entropy at a fixed fraction of the virial radius and the temperature, such 
that S oc T ' 65 (see the left-hand panel of Fig. 0J). This slope is in agreement with that expected 
from the observed scaling of X-ray luminosity with temperature. The equivalent entropy-mass 
relation is S oc M 0,36 , again shallower than expectations from simple gravitational collapse 
( Pratt et allbOOfih . 




Figure 4: Left: Entropy measured at 0.1 -R200 vs system temperature for 66 systems ffrom lPonman et al.ll2003l) . 
The dotted line shows the relation expected from pure gravitational structure formation (S oc T). The solid line 
is the best-fitting power law relation: S oc y - 65 . Right: Entropy profiles of 10 c lusters spanning t he mass range 
1 - 12 x lO 14 /!™ 1 Mq (~ 2 [blue] - 8 keV [red]), observed with XMM-Newton Cfrom lPratt et allkOOrJ) . The profiles 
have been scaled using the best-fitting empirical entropy scaling. 



High resolution entropy profiles can give essential insights into the physical processes at 
play. The right-hand panel of Fig. 0] shows the entropy profiles of the same ten clusters from 
Fig |31 this time res caled by the observed empirically-determined entropy-temperature relation 
(S oc T 65 ; see also Piffaretti et al. . 200.4 ) . Two important points are evident from this exercise. 
Firstly, outside the core regions, the entropy profiles show a remarkable degree of similarity, 
and the radial dependence of entropy is well described by the relation S°$R , as expected 
from standard gravitational collapse (Sect. I^J), and as is seen in gravitation-only simulations 



( Voit et all 120 05). Thus non-gravitational processes appear to alter the power-law scaling of the 



entropy-temperature relation without breaking structural similarity, at least for clusters of this 



mass range and above (results for the group regime < 2 keV are less conclusive, see lMahdavi et al 



2005). It is still unclear what mechanism is causing this effect. One proposal is that smoothing 
of the gas density due to preheati ng in filaments and /or infalling group s may boost entropy 
production at the accretion shock ( Voit et al. . 200.4 Ponman et all 2003). However, numerical 
simulations run specifically to test this effect seem to suggest that the boo sting, though present , 
appears to be substantially diminished in the presence of radiative cooling (Borgani et al ., 2005). 
The same simulations suggest that feedback from supernovae is too localised to have a significant 
effect on smoothing of the accreting gas. Another possibility is that activity from the central 
AGN produces the extra heating, although the observed normalisations require than the AGN 
affects the entropy distribution in a non-catastrophic manner at least out to i?iooo(~ 0.5i?20o) 
and it is not clear why and how this would produce the observed power-law entropy-temperature 
dependence. 



Secondly, inside the core regions the dispersion increases dramatically, constituting an ef- 
fective breaking of similarity. Three probable mechanisms for breaking of similarity in the core 
are radiative cooling, AGN feedback, and mixing of high and low entropy gas due to merging 
activity. Disentangling the effects of these processes will require a concerted effort from both 
observations and simulations. 



6 Summary and perspectives 

The explosion of high quality data available from the most recent generation of X-ray satellites 
has allowed unprecedented insights to be gained into cluster physics. The parallel increase 
in computing power available for large hydrodynamic simulations has facilitated increasingly 
sophisticated comparisons between observations and theory. The resulting cross-fertilisation 
has had profoundly beneficial effects on our understanding of the formation and evolution of 
structure in the Universe. 

X-ray observations have allowed strong constraints to be put on the form of the mass density 
profile in clusters, suggesting a cusped form and a variation with mass in agreement with the 
predictions from numerical simulations. This result is true over the entire measurable radial 
range (~ 0.001-0.8 -R200) m morphologically relaxed clusters, groups, and individual galaxies. 

The observed X-ray mass-temperature (M-T) relation has been measured using data from 
XMM-Newton and Chandra with several samples of moderate size, yielding results which are in 
very good agreement both in terms of slope and normalisation. This agreement can be attributed 
both to better data quality and to improved data analysis techniques. However, there is still an 
offset between the observed M-T normalisation and that determined from numerical simulations. 
While a part of this discrepancy may come from details of the X-ray analysis (perhaps principally 
the assumption of hydrostatic equilibrium), it is also clear that since simulations currently do 
not reproduce the gross scaling properties of the cluster population (e.g., the Lx~T relation) 
the physics governing the baryonic component are incompletely understood. 

We can get an insight into this by looking at the entropy of the ICM. The entropy-temperature 
relation is shallower than expected, S oc T - 65 , such that poor systems have higher entropy than 
rich systems relative to that expected from pure gravitational collapse. Entropy profiles of clus- 
ters formed purely through gravitational processes are self-similar (once scaled appropriately) 
and are characterised by a radial dependence S oc R 11 outside the central regions (R > 0.1 -R200)) 
a slope characteristic of shock heating of the ICM. The entropy profiles of observed clusters are 
self-similar outside the central regions once scaled by the empirical scaling, and display a slope 
similar to that expected from shock heating. Thus the entropy is higher throughout the ICM in 
poor systems. It is not yet clear what physical mechanism is responsible, but the most obvious 
candidates are preheating and AGN activity. In the central regions (R < 0.1i?20o)> there is a 
clear breaking of similarity. This is probably due to the combined action of radiative cooling, 
AGN heating and mixing of high and low entropy gas during mergers. 

Most of the above results have been derived from observations of morphologically relaxed 
systems, and are thus they are neither representative of the cluster population as a whole, nor 
useful from the point of view of establishing the intrinsic scatter about the scaling relations. 
The challenge now is in the observation of larger samples with unbiased mass and redshift 
sampling. These should be used to i) further test the self-similarity of form; ii) derive the 
exact slope, normalisation and intrinsic scatter of the M-T relation, requiring cross-calibration 
and combination of different mass estimation methods including X-ray, lensing, and dynamical; 
iii) compare with state of the art simulations to probe the sources and timescales of the non- 
gravitational energy input which alters the similarity of the cluster population. 
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